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Using the known solutions of the Lugiato-Lefever equation, we derive universal trends of Kerr frequency combs. In particular, 
normalized properties of temporal cavity soliton solutions lead us to a simple analytic estimate of the maximum attainable 
bandwidth for given pump-resonator parameters. The result is validated via comparison with past experiments encompassing a 
diverse range of resonator configurations and parameters. © 2013 Optical Society of America 
OCIS codes: 230.5750, 190.5530, 190.4380 



Over the last few years, the generation of frequency combs 
in high-Q Kerr microresonators pumped by continuous-wave 
(cw) laser light has attracted considerable interest JT]. De- 
spite this interest, the dependence of comb characteristics on 
the pump-resonator parameter-space remains largely unex- 
plored, with little or no consensus existing on guidelines for 
Kerr comb optimization. This shortcoming originates from 
the computational complexity of traditional models, such as 
the coupled-mode equation model JU^J. Approximate ana- 
lytic solutions exist but the comb characteristics still cannot 
be inferred in closed form . 

We have recently shown that Kerr frequency combs can be 
efficiently modeled using a generalized mean-field Lugiato- 
Lefever equation (LLE), and that they can be associated with 
the cavity soliton (CS) solutions of this equation (6). Also 
Herr et al. have presented convincing experimental evidence 
in strong support of the CS hypothesis Q. Here, we use the 
mean-field framework to identify universal trends in the dy- 
namics and characteristics of Kerr frequency combs. Our ob- 
jectives are two-fold: (i) we link the known solutions of the 
LLE into the Kerr comb context and (ii) use the solutions to 
derive universal scaling laws that allow comb bandwidths to 
be analytically estimated. Comparison with past experiments 
across a wide variety of resonator parameters and architec- 
tures shows good agreement with our results. 

Our starting point is the normalized LLE ||6]|8], 
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Here t is the slow time describing the evolution of the intra- 
cavity electric field envelope E(t,x) at the scale of the cav- 
ity photon lifetime while T is a fast time defined in a refer- 
ence frame traveling at the group velocity of light in the res- 
onator. The terms on the right-hand side of Eq. ([T) describe, 
respectively, cavity losses, Kerr nonlinearity, pump-cavity de- 
tuning, 2nd-order group-velocity dispersion (GVD), and ex- 
ternal pumping. 77 is the sign of the GVD coefficient fc, and 
we assume here anomalous dispersion, 77 = — 1 . The normal- 
ization follows Ref. (8): t -> at/t R , x -> T(2a/|/3 2 |L) 1 / 2 , and 
E -> E{yL/a) l l 2 . Here f R = FSR 1 is the cavity roundtrip 
time with FSR the free-spectral range, a = (a\L+ 9)/2 is 
the total round-trip loss parameter with a; describing internal 
linear absoiption and 9 the input coupler power transmission 



coefficient, and L is the resonator length. The nonlinearity co- 
efficient 7 = 27T«2/(ApA e ff) with «? the nonlinear refractive 
index, X p the pump wavelength, and A er r the effective mode 
area. The solutions of Eq. ([TJ are governed by two normalized 
parameters: the pump strength S = E m (jL0 / a 3 ) 1 ? 2 , with E m 
the cw pump amplitude in units of W 1 / 2 , and the detuning 
A = 5()/a, where 5q is the phase detuning of the pump from 
the closest cavity resonance. Note that A is formally identical 
to the £0 parameter introduced in Ref. . 

Under cw pumping, the simplest steady-state solutions of 
Eq. (fl} are homogeneous (dE/dx = 0) and satisfy 



X = Y 3 - 2AY 2 



(A 2 +1)T, 



(2) 



where X = \S\ 2 and Y = \E\ 2 are the normalized pump and in- 
tracavity powers, respectively. Equation d2} is the well-known 
cubic equation of dispersive optical bistability J5). At con- 
stant pump, it also describes the Kerr tilt of the cavity reso- 
nances as shown in Fig. 02a) for X = 10 (black curve). Here 
X is large enough for Y(A) to be multivalued. The lower 
branch, existing for A > A T , is always stable [A-j- ~ 3(X/4) 1 ^ 
for X> 1] while the intermediate branch (dotted black) is 
homogeneously unstable. The left side of the resonance, in- 
cluding the upper branch of the multivalued region, exhibits 
modulation instability (MI) for intracavity powers above the 
threshold Y > 1 [9], or equivalently for detunings A > Ami = 
1 — \/X — 1 (dashed black). Accordingly X — 1 is the absolute 
MI threshold, also referred to as the hyperparametric thresh- 
old in the comb literature 0~|-[3] . 

MI leads to patterned solutions J9) and the green curve in 
Fig. [TJa) represents the peak power of a range of such solu- 
tions (X = 10 as above). They are steady-state periodic solu- 
tions of Eq. ([T]) and have been obtained with a Newton solver 
and a continuation method J6| . Here the MI frequency is set to 
the most unstable one at threshold. For the point highlighted 
by the green cross, we have also plotted in Figs. Qlb)-(c) the 
corresponding temporal and spectral intensity profiles. As can 
be seen, an MI solution corresponds in the spectral domain to 
a frequency comb, but resonators are typically long enough to 
fit multiple periods of the pattern, so the comb modes would 
be separated by multiple FSRs. For larger A, the MI branch 
becomes unstable (dotted green), and split-step simulations 
of Eq. (Q3 reveal a transition to chaotic combs whose broad- 
linewidth components are separated by a single FSR, as also 
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Fig. 1. (Color online) (a) Intracavity peak power versus detuning 
for the cw (black), MI (green), and CS (blue) branches. Dashed and 
dotted lines indicate unstable solutions as described in the text. The 
temporal and spectral intensity profiles of particular MI and CS so- 
lutions (crosses) are illustrated in (b,c) and (d,e), respectively. 

observed from the coupled-mode equation model |0J. 

CSs constitute the last class of steady-state solutions of 
Eq. (Q~|) J6][8l. Their peak power, obtained as for the MI 
branch, completes Fig. [Ha) [blue curve]. In the time domain, 
CSs are localized pulses sitting atop a cw background and, 
due to resonator periodic boundary conditions, correspond to 
a frequency comb in the spectral domain J6]. A typical so- 
lution (blue cross) is illustrated in Figs. [T{d)-(e). Here the 
comb separation is a single FSR as there is only one CS in the 
cavity. The CS background matches with the lower state ho- 
mogenous solution, as it is the only stable cw solution, hence 
CSs only exist above the up-switching point, A > A+. It is also 
found that CS existence is limited to a maximum detuning of 
about n 2 X /8 ifTol . The lower part of the CS branch is always 
unstable (dotted blue). On the upper branch, CSs go through a 
Hopf bifurcation for lower detunings (dashed blue). Here CSs 
are oscillating (breathing) over multiple round-trips [H0UT2I . 
More complex chaotic regimes also exist in this region |[T2"1 . 

In the context of Kerr frequency combs, important conclu- 
sions emerge from the (co)existence of the cw, MI, and CS 
solutions represented in Fig.Q~}a). CSs are the preferred solu- 
tions for combs because, in comparison to MI, they are stable 
over a wider parameter range and their peak power is higher, 
leading to broader spectral bandwidths. However, CSs can- 
not spontaneously emerge from an initial cw background as 
stable CS solutions are disconnected from the cw branch fl8]. 
Moreover, to exploit thermal self-locking, the resonance is 
approached from the high frequency side [from the left in 
Fig.Q2 a )]- Combined with the fact that Ami is always smaller 
than A-s, it implies that MI will always occur first. The onset of 
MI can therefore be associated with the "primary combs" de- 
scribed in other works [2 |. The MI branch extends into the CS 
branch, but unstable (chaotic or otherwise) MI and CS states 
lie in between the onset of MI and the desirable stable CSs. 
Although beyond the scope of this Letter, preliminary sim- 
ulations show that a chaotic MI state can condensate into a 
set of CSs with appropriate ramping of A. Note that the chain 
stable MI — > chaotic MI — > stable CSs as implied by Fig.Qla) 
is precisely that observed in recent experiments Q. 

In the CS regime, the comb bandwidth is determined by 
the temporal duration of the CS. To gain some general in- 
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Fig. 2. (Color online) (a) Normalized CS duration (left axis) and 
peak power (right axis) versus A/X for various pump levels (dashed: 
oscillatory CSs; other unstable CSs are omitted for clarity), (b) Ex- 
amples of numerical CS temporal intensity profiles for parameters 
corresponding to the dots in (a). In both panels, the numerical re- 
sults are compared with analytic approximations (grey circles). 

sights, we have calculated CS characteristics across a wide 
range of parameters. In Fig. |2ja), we plot the full-width at 
half-maximum (FWHM) At of the CS solution (left axis) and 
its peak power Fpeak (right axis) as a function of cavity de- 
tuning for several values of X. Note how we have normal- 
ized all axes based on X to reveal universal CS trends. We 
see that the CS peak power increases linearly with detun- 
ing, with a slope approaching two for increasing values of 
X, i.e., Fpeak — 2A. As the peak power increases, the CS dura- 
tion decreases and approaches At ~ 1 .763 /VA for X 1 . In 
Fig.|2b) we also show normalized temporal intensity profiles 
of selected solutions [dots in Fig. 12a)] superimposed with a 
sech 2 pulse with peak power and duration as derived above. 
Agreement is excellent for all parameters, implying that CSs 
can be approximated as £cs — v / 2Asech(v / AT). It is worth 
noting that this expression is an exact analytical solution of 
the normalized LLE (Q]i for a pulsed pump identical to the 
soliton, 5(t) = £cs, and is also found by perturbation the- 
ory as the fixed point of (fl} JTE], as well as in the solution 
derived by Barashenkov and Smirnov with X = A ^> 1 I5lfl0ll . 

Considering £cs to be a good approximation of the 
intracavity field allows us to derive a simple estimate for 
the comb bandwidth obtainable for given pump-resonator 
parameters. Specifically, assuming on-resonance pumping 
(A = X) and critical coupling (a = 9), and transforming 
into physical units, we get the following simple theoretical 
estimate for the 3-dB comb bandwidth, 

_0J15 / 2yP m eApFSR _ 0.315 \ 2yP m & 
/theo " 1.763 V rcc|j3 2 | ~ 1.763 V 7r|j3 2 | ' W 

with P; n = |£in| 2 the pump power, Q the quality factor of the 
cavity, and c the speed of light. Eq. (01 also reveals the rel- 
evance of the various parameters. Recalling that the cavity 
finesse & ' = n/a = QXpFSR/c, it is clear that the suitability 
of a given resonator for broadband comb generation is de- 
termined solely by its finesse (or losses), nonlinearity, and 
GVD coefficient, and is independent of its FSR. Note that MI 
generally follows trends similar to CSs such that Eq. (|3) is 
qualitatively correct irrespective of the operation regime. 
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Table 1. Comparison of experimental 3-dB comb bandwidths (A/ eX p) with analytical predictions (A/ t h eo )- 



Ref. 


/[W-'km- 1 


Pin [mW] 


Xp [nm] 


Q 


FSR [GHz] 




ftlps^km- 1 


A/exp [THz] 


A/ th eo [THz] 





0.405 


24 


1553 


400- 10" 


35.2 


73 • 10 J 


-5.9 


1.6 


1.59 


m 


15 


160 


1550 


100 ■ 10 6 


882 


460- 10 3 


-6.3 


31 


36.4 


ma 


800 


2000 


1562 


3 ■ 10 5 


226 


350 


-47 


10 


15.6 


na 


0.032 


55.6 


1560 


1.90 -10 9 


18.2 


180- 10 3 


-13 


0.6 


0.71 


El 


25 


1000 


1550 


270- 10 6 


850 


1.2 • 10 6 


-4.0 


40 


388 



The validity of Eq. 01 has been checked by comparing its 
predictions against experimental comb bandwidths for a num- 
ber of references, as summarized in Tabled] with all relevant 
experimental parameters also listed. Aside from the comb re- 
ported in ifTTl . agreement is good throughout, with discrep- 
ancies attributed to inaccuracies in experimental parameters, 
difficulties in estimating the bandwidths of highly-structured 
experimental combs, as well as approximations implicit in 
Eq. (0. Note in particular how the comparatively low finesse 
and high GVD of the resonator used in 0151 is offset by the 
large nonlinearity {yP m product). The large deviation between 
the estimated and experimental bandwidths for IfTTl stems 
from the LLE (0 and Eq. (0 not taking into account higher- 
order dispersion. This approximation becomes inaccurate for 
comb spectra extending into the normal GVD regime, and re- 
sults in significant over-estimation of the bandwidth. We can 
mitigate this issue by recalling that CSs perturbed by higher- 
order dispersion emit dispersive waves (DW) into the normal 
GVD regime ||6]|T8][19). Assuming that the 3rd-order GVD 
coefficient fa is the dominant contribution, the DW frequency 
shift can be written as A/bw = 3 1 JS2 1 / (2 tt/33 ) . Because the 
generated DWs are linear (they reside in a spectral region 
not supporting solitons), no significant spectral broadening 
is expected beyond their spectral shift. We can thus interpret 
A /iim = 2A/dw as being an upper limit of the full- width band- 
width of a Kerr comb. Interestingly, while Eq. (0 suggests 
that low GVD is beneficial for broadband combs, the full- 
width limit bandwidth Afn m is in this case constrained to a 
lower value. It is only a low dispersion slope ^3 (i.e., flat dis- 
persion) that can lift this limit. These trends are illustrated in 
Fig. where we plot the bandwidths estimated from Eq. (0 
as a function of fa for various values of the product yP m ^ 
(solid lines) as well as the full-width upper limit A/i; ra (dot- 
ted line), assuming here fa =0.1 ps 3 /km. In Fig. we also 
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Fig. 3. (Color online) Theoretical Kerr comb bandwidth calculated 
from Eq. ((5} versus 2nd-order GVD for varying values of jP m ,^ 
(solid lines) and upper limit (dotted line). The black diamonds are 
the theoretical estimates for the References listed in TablefTJwith the 
corresponding experimental bandwidths shown as green dots. 



show the experimental (green dots) and theoretical [Eq. (0; 
black diamonds] bandwidths extracted from Table Note in 
particular how the comb bandwidth observed in JT7], overes- 
timated by Eq. (0, falls on the theoretical upper limit. 

To conclude, we have analyzed the solutions of a dimen- 
sionless mean-field equation from the Kerr comb perspec- 
tive. By identifying universal trends in the solutions we have 
derived universal scaling laws linking experimental pump- 
resonator parameters to the obtainable comb bandwidth. 

We acknowledge support from the Marsden Fund (Royal 
Society of New Zealand) and useful discussions with T. Herr. 
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